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ABSTRACT 

We use N-body simulations to study the tidal evolution of globular clusters (GCs) 
in dwarf spheroidal (dSph) galaxies. Our models adopt a cosmologically motivated 
scenario in which the dSph is approximated by a static NFW halo with a triaxial 
shape. We apply our models to five GCs spanning three orders of magnitude in stellar 
density and two in mass, chosen to represent the properties exhibited by the five 
GCs of the Fornax dSph. We show that only the object representing Fornax's least 
dense GC (Fl) can be fully disrupted by Fornax's internal tidal field — the four denser 
clusters survive even if their orbits decay to the centre of Fornax. For a large set of 
orbits and projection angles we examine the spatial and velocity distribution of stellar 
debris deposited during the complete disruption of an Fl-like GC. Our simulations 
show that such debris appears as shells, isolated clumps and elongated over-densities 
at low surface brightness (> 26 mag/arcsec 2 ), reminiscent of substructure observed 
in several MW dSphs. Such features arise from the triaxiality of the galaxy potential 
and do not dissolve in time. The kinematics of the debris depends strongly on the 
progenitor's orbit. Debris associated with box and resonant orbits does not display 
stream motions and may appear "colder" / "hotter" than the dSph's field population if 
the viewing angle is perpendicular/parallel to progenitor's orbital plane. In contrast, 
debris associated with loop orbits shows a rotational velocity that may be detectable 
out to few kpc from the galaxy centre. Chemical tagging that can distinguish GC debris 
from field stars may reveal whether the merger of GCs contributed to the formation 
of multiple stellar components observed in dSphs. 

Key words: galaxies: halos - Galaxy: evolution - Galaxy: formation - Galaxy: 
kinematics and dynamics 



1 INTRODUCTION 

The sizes (rhalf ~ 100 pc), velocity dispersions (cry ~ 10 km 
s _1 ) and luminosities (Ly ~ 10 3_7 Z/©) of dwarf spheroidal 
(dSph) galaxies imply that these objects are the most dark- 
matter dominated galaxies known (see reviews by Mateo 
1998 and Gilmore et al. 2007). Their dynamical mass-to- 
light ratios are sufficiently large (M/L v ~ 1O 1_3 [M/Ly] ) 
that their stars can be treated as test particles tracing their 
underlying gravitational potentials. Dynamical constraints 
on the dark matter distribution in dSphs may therefore be 
compared directly with the predictions from cosmological 
simulations to study the elusive nature of dark matter (e.g. 
Gilmore et al. 2007 and references therein). 

In the last few years, however, wide-field imaging and 
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multi-object spectroscopic surveys have revealed that dSph 
stellar components are more complex than previously be- 
lieved. Localized oddities such as shells (Coleman et al. 2004, 
2005; Olszewski et al. 2006), butterfly-shaped substructures 
(Stetson & Hesser 1998) and cold clumps (Kleyna et al. 
2003, 2004; Ibata et al. 2006) have been reported in sev- 
eral dSphs. These results conflict with the naive expectation 
that dSph dark matter haloes (with typical crossing times of 
a few hundred Myr) should host well-mixed stellar popula- 
tions and may pose a challenge to the present cosmological 
paradigm: Cold Dark Matter (CDM). In a CDM cosmogony, 
dark matter haloes follow a density profile that is nearly 
independent of their mass (Navarro, Frenk & White 1996, 
1997; hereafter NFW) and scales as p oc r~ 7 in the inner- 
most regions. Navarro et al. (2008) have recently shown that 
7 varies slightly with halo mass and that dark matter halo 
profiles only deviate slightly from NFW models (7 = 1). 
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In such "cuspy" potentials, unbound stellar substructures 
would dissolve after a few crossing times, as shown by the 
N-body simulations of Kleyna et al. (2003), who also show 
that the lifetimes of these substructures would be consid- 
erably prolonged if the halo has a constant-density "core". 
Thus it has been argued that the presence of stellar sub- 
structures in dSphs is in apparent contradiction with CDM 
expectations. 

A second oddity in this context is the existence of glob- 
ular clusters (GCs) in dSphs. Although GCs are consider- 
ably more common in dwarf irregular galaxies (e.g. Georgiev 
ct al. 2008), at least two of the Milky Way dSphs — Fornax 
and Sagittarius — contain GCs (five and four, respectively). 
Among the Milky Way dSphs these galaxies are also the 
largest and most luminous (Mateo 1998), which may be re- 
lated to their ability to form GCs at an early epoch. The 
extended spatial distribution of Fornax's GCs (most lie, in 
projection, near Fornax's half-light radius of ~ 400 pc) has 
fueled arguments against cuspy NFW profiles in dSphs. For 
example, Goerdt et al. (2006) demonstrate that dynamical 
friction, proportional in strength to the background density, 
requires only a few Gyr to drag GCs from Fornax's half-light 
radius to its centre in a cuspy potential. Oh, Lin & Richer 
(2000) note a possible solution; namely, that strong tidal in- 
teractions between the Milky Way and Fornax could inject 
energy into the GC orbits and delay their decay. However, 
this scenario is now ruled out by proper motion measure- 
ments (Dinescu et al. 2004; Piatek et al. 2007; Walker et 
al. 2008) that imply Fornax follows a nearly circular or- 
bit around the Milky Way. Hernandez & Gilmore (1998), 
and also more recently Goerdt et al. (2006) and Read et al. 
(2006), pose a second viable solution, showing that dynami- 
cal friction is suppressed in the inner regions of cored halos, 
thereby providing a natural explanation for the apparent 
levity of Fornax's GCs. The survival of GCs in dwarf galax- 
ies may then provide important clues regarding the inner 
structure of dark matter haloes. 

Despite recent theoretical efforts to understand the for- 
mation mechanism of GCs in a cosmological context (e.g. 
Prieto & Gnedin 2008, Bournaud et al. 2008), the forma- 
tion of GCs in dwarf spheroidals lacks definite theoretical 
underpinning. For example, it is unclear whether Fornax 
and Sagittarius are the only MW dSphs that satisfy the re- 
quired conditions to host GCs or whether GCs also formed in 
smaller and fainter dSphs but did not survive to the present 
epoch. The present contribution elaborates on the latter sce- 
nario by studying the survival of Globular Clusters in the 
tidal field of a typical dwarf galaxy and attempts to address 
the signatures that the tidal disruption of stellar clusters 
may have left in the stellar populations of dSphs. We ex- 
amine these issues here with the aid of N-body simulations 
that adopt a CDM-motivated dSph model, where the dwarf 
galaxy potential is assumed to follow an NFW profile with 
triaxial shape. This is a timely issue in light of ongoing spec- 
troscopic surveys involving thousands of target stars in Local 
Group dSphs. 

This paper is organized as follows. Section 2 introduces 
the observational data and the modelling procedure. Sec- 
tion 3 describes the numerical technique and its limitations. 
In Section 4 we study the mass and orbital evolution of 
GCs in NFW halos. In Section 5 we examine the spatial 
and velocity distribution of disrupted GC debris. Section 6 



discusses the possible detection of debris in dSphs. We end 
with a brief summary in Section 7. 



2 MODELS AND OBSERVATIONAL DATA 

Our models assume that the mass profile of a dSph can be 
approximated by a triaxial Navarro, Frenk & White (1996, 
1997) dark matter halo. Due to the high mass-to-light ra- 
tios estimated in dSphs (e.g. Mateo 1998, Simon & Geha 
2007, Penarrubia, McConnachie & Navarro 2008) stars are 
assumed to contribute negligibly to the overall potential. 

Our GC models are N-body realizations of a King 
(1966) profile orbiting within a static NFW potential. The 
masses of GCs are relatively small in comparison with the 
host halo mass and here we assume that the presence of glob- 
ular clusters (GCs) does not alter the distribution of dark 
matter. Dynamical friction resulting from two-body inter- 
actions between the GC and the background dark matter 
particles is introduced in our computations as an external 
friction term following the method outlined in Penarrubia 
et al. (2006) (see §[3]). We also assume that the dSph models 
are in isolation, thereby neglecting any effects that external 
forces (e.g. the Milky Way tidal field) may induce on the 
orbital evolution of the GC. 

For illustrative purposes, our study focuses on the dy- 
namical evolution of a GC population chosen to resemble 
that of the Fornax dSph. Fornax hosts the largest number 
of known GCs (five) among all Local Group dSphs. Insofar 
as the scatter in parameters characterizing the dark matter 
halos surrounding MW dSphs is small (e.g., Penarrubia et 
al. 2008a, b), and given the wide range of masses and spa- 
tial sizes covered by Fornax's cluster population, our results 
should capture the most relevant features of the evolution 
of GCs in dwarf galaxies. In the remainder of this Section 
we provide further details on the numerical methods. 



2.1 The dwarf galaxy model 

We assume that the dSph galaxy is embedded in a static 
dark matter halo that follows an NFW density profile, 

Po 



Pnfw 



{r/r s )(l + r/r s ) 2: 



(1) 



Po = 



V, 



max' max 



47rGri[ln(l + r max /r s ) - r mB x/r s /(l + r ma x/r s )] ' 



where r 3 the scale radius and Knax = Vc(r max ) is the peak 
circular velocity. For NFW halos, r max ~ 2.17r s . The halo 
mass at a given radius is simply M(r) — 47rrfpo[hi(l + 
r/r 3 ) - r/r„/(l + r/r s )]. 

In the Cold Dark Matter paradigm, halos generally have 
a triaxial shape (Davis et al. 1985, Frenk et al. 1998, Jing 
et al. 1995, Jing & Suto 2002) that depends only weakly on 
radius (Kazantzidis et al. 2004b, Allgood et al. 2006). We 
account for that effect by introducing elliptical coordinates, 
using the transformation r — > m, where 



2 2 2 

X X z 

a 2 b 2 c 2 ' 



(2) 



and (a, b, c) are dimensionless quantities. Our coordinate 
frame is selected so that a > b > c. We have used the formu- 
lae of Chandrasekhar (1969) to solve the Poisson equations 
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in elliptical coordinates for the density profile of eq.|2|. The 
resulting potential can be written as 



$nfw(i, y, z) = 27vGabcp r s (3) 
TriM/r. dr 

+ const, 



projected surface density 



1 + m(r)/r s Va 2 + ">V& 2 + ">Vc 2 + • 



where m 2 (r) = x 2 /(a 2 + r) + y 2 /{b 2 + r) + z 2 /(c 2 + r). 
We choose a triaxiality of (a, 6, c) = (1.47, 1.22, 0.74), which 
roughly corresponds to the most common shape found in 
cosmological simulations (e.g. Jing & Suto 2002). The re- 
cent work of Allgood et al. (2006) shows that haloes tend 
to become slightly rounder as their mass decreases and also 
at outer radii. Because GCs in dSphs are only found in the 
inner-most regions of the dwarf dark matter haloes (see be- 
low), the triaxiality gradients found by Allgood et al. (2006) 
can be safely neglected for this study. Note also that choos- 
ing a different combination of axis-ratios will change the 
phase-space volume occupied by each orbital family (see e.g. 
Fig. 1 of Adams et al. 2007) without affecting the qualita- 
tive analysis on the dynamical evolution of GCs in a triaxial 
potential outlined in the following Sections. 

The two free parameters of the NFW potential (V max 
and r max ) are scaled to the values estimated for MW dSphs 
by Penarrubia et al. (2008a,b) who find a remarkably small 
scatter in dSph masses. Indeed, the studied dSphs differ by 
less than a factor ~ 4 in V max , with peak velocities in the 
range 10 < Vmax < 40 km/s and peak radii 2^-r max ^9 
kpc, despite the fact that their luminous components span 
four decades in luminosity (~ 10 3 -10 7 Lq) and one decade in 
spatial size (100-1500 pc core radii). Systems of particular 
interest for this contribution are Sagittarius and Fornax, as 
these are the only MW dSphs with associated populations of 
gravitationally bound GCs. The estimated NFW parameters 
of these galaxies are, respectively, Knax = 13.8 and 20.6 
km/s, with peak velocity radii at r max = 2.6 kpc and 4.1 
kpc. 

Finally, we note that in our dSph models stars and GCs 
are confined to the inner-most regions of the halo (i.e. the 
"cusp" in an NFW halo, see Penarrubia et al. 2008a). Due 
to (i) the small amount of dark matter substructures in sub- 
haloes with respect to that in the main halo and (ii) the large 
galactocentric distances at which substructures tend to lo- 
cate (Springel et al. 2008), for simplicity we neglect here 
the effects that the halo granity may induce on the orbits 
of cluster particles and assume a smooth potential for our 
dwarf halo models. 



T,(R) = k{ 



[l + {R/Rc) 2 ] 1/2 [l + (Rt/Rc 



1/2 



(5) 



where R is the projected radius, k is related to the central 
surface brightness by Eo = fc^l — 1/^/1 + cj^j and ck = 
Rt I Rc ■ 

Table [l] lists the properties of the GC population of the 
Fornax dwarf relevant for this study. For comparison we also 
include the properties of the four clusters in the Sagittarius 
dwarf. Throughout this work we convert luminosities into 
masses and vice versa assuming GCs have stellar mass-to- 
light ratios of approximately (M/L)* w 3 (e.g. Bruzual & 
Chariot 1993, Bell & de Jong 2001). As we show below, the 
orbital and dynamical evolution of a GC due to dynamical 
friction and tidal disruption depend primarily on the GCs 
mass and density. 

Fig. [T] compares the density profiles of the Fornax GCs 
to those of the NFW halos for the whole range of Vm ax val- 
ues estimated for the Milky Way dSphs by Penarrubia et 
al. (2008a). In particular, filled circles show the density pro- 
file of the Fornax halo. This figure demonstrates two rele- 
vant points. First, Fornax's GCs divide into three categories 
based on their central densities: (i) low-density clusters like 
Fl, with p*(0) ;$ Pnfw(-Rc); (ii) intermediate-density clus- 
ters (F2), p*(0) <; Pnfw(-Rc) and (iii) high-density clusters 
like F3, F4 and F5, with p*(0) S> Pnfw(Rc)- Henceforth, to 
simplify our study we will consider only the clusters Fl, F2 
and F3 to represent the three GC categories. Second, the 
NFW halos estimated for the MW dSphs display a small 
scatter in density at small radii (r C r,). Varying the peak 
velocity between 10 and 40 km/s leads to a density variation 
of a factor £3 at a given radius (e.g. at 100 pc), whereas 
the stellar densities of the Fornax and Sagittarius GCs span 
a factor ~ 1000. Then the domination of dSph gravitational 
potentials by similar dark matter halos implies that the evo- 
lution of GCs depends primarily on the clusters' properties 
at the time of formation, rather than on the properties of 
their host halos. 

The second parameter relevant to this study is the mass 
of the GC. Taking into account that all the bright dwarfs 
in the Milky Way appear to have dynamical masses close 
to M dyn ~ 1O 7 M (e.g. Mateo 1998), GCs represent only 
a small fraction of the dwarf mass, with mass ratios in the 
range Mac/M dyn ~ 10 -3 (Fl) and 10~ 2 (F3). 



2.2 Globular Cluster models 

Our N-body realizations of GCs assume a King (1962) stellar 
density distribution, 



P* = 



K 



1 + (r/r e ) 2 
1 + (rt/r c ) 2 



1/2 



(4) 



where r c and rt are the King and tidal radii, and K is an 
arbitrary constant that we select as K = L/ J Q r * p*d 3 r, where 
L is the total luminosity of the cluster. All our cluster models 
have a particle number N = 10 5 . 

For comparison with observations it is easier to use the 



3 NUMERICAL SETUP 
3.1 The N-body code 

We follow the evolution of the cluster N-body models in the 
dwarf potential using Superbox, a highly efficient particle- 
mesh gravity code (see Fellhauer et al. 2000 for details). 
Superbox uses a combination of different spatial grids in 
order to enhance the numerical resolution of the calculation 
in the regions of interest. In our case, Superbox uses three 
nested grid zones centered on the highest-density particle 
cell of the cluster. This centre is updated at every time step, 
so that all grids follow the cluster along its orbit. 

Each grid has 64 3 cubic cells: (i) the inner grid has a 
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log 1D r' (kpc) 

Figure 1. Density profiles of the five GCs in Fornax (F1,...,F5) 
and of the underlying dark matter halo. The halo of the For- 
nax dSph is assumed to follow an NFW profile with Vmax — 20 
km/s. For comparison we also plot the density profiles of halos 
with V m ax = 10 km/s and 40 km/s, which represent the approx- 
imate range of maximum circular velocities of the Local Group 
dSphs (Penarrubia et al. 2008b). We distinguish three categories 
of clusters: Those with central densities higher (F3, F4 and F5), 
comparable (F2) and lower (Fl) than that of the dark matter 
halo at r ~ R c . 

spacing of dx± = 10i? c /62 ~ 0.16ii c and is meant to re- 
solve the entire Globular Cluster, (ii) The middle grid ex- 
tends to intermediate distances with a resolution of dx2 = 
100-Rc/62 ~ l.GRc- (iii) The outermost grid extends out to 
250 kpc and is meant to follow particles that are stripped 
from the cluster and subsequently orbit within the main 
galaxy. 

In addition to the overall dSph potential, GCs feel dy- 
namical friction as they orbit within the dSph's extended 
halo. To describe dynamical friction in aspherical haloes 
we apply the equations of Binney (1977), which recover 
the Chandrasekahr (1943) formula for a = b = c = 1. 
Penarrubia, Just & Kroupa (2004) test the accuracy of 
these equations against self-consistent N-body simulations 
and find that, for a large range of orbital parameters and 
masses, the overall radial evolution of point-masses within 
flattened dark matter haloes can be reproduced to within 
a few percent. As the N-body simulations of Penarrubia et 
al. (2002) show, halo asphericity introduces a strong depen- 
dence between the orbital decay time and the orientation 
of the orbit with respect to the halo's plane of symmetry. 
In particular, these authors show that orbits that lie on the 
major-intermediate axis plane have shorter decay times. 

For simplicity we assume here that the velocity distri- 
bution of the dark matter halo is Maxwellian and isotropic. 
Kazantzidis et al. (2004a) have examined the inadequacies of 
the local Maxwellian approximation when applied to NFW 
haloes and, although it may introduce problems to gener- 
ate haloes in perfect equilibrium, the deviations from Gaus- 



sianity of the exact distribution function have a negligible 
impact on the friction term (see their Fig. 4). Under these 
approximations, the friction force per unit mass can be writ- 
ten as (Binney 1977) 

f d( = -2V2^G* / 9nfw^BM gc In A, (6) 

where cr 2 (m) = /0^ FW J*°° dr' M{r')pnwr{i'') /r 2 is the halo 
velocity dispersion, B — dr exp[— v 2 /(l + r)/2a 2 ]/(l + 
t) 5,/2 and Mgc is the cluster's bound mass. Penarrubia et al. 
(2004) find that the best fit to self-consistent N-body orbits 
is obtained for a Coulomb logarithm In A = 2.1, in good 
agreement with the recent results of Arena & Bertin (2007). 
However, it is well known that this quantity is fairly sensitive 
to numerical aspects, such as particle number and spatial 
resolution (e.g. Prugniel & Combes 1992, Wahde & Donner 
1996) as well as to the host's density profile (Hashimoto et 
al. 2003, Just & Penarrubia 2005, Read et al. 2006, Arena 
& Bertin 2007). It is therefore necessary to note that the 
decay times of clusters orbiting in dark matter haloes are 
uncertain to a significant degree, since orbital decay rates 
are proportional to the value assumed for In A. 

In order to implement eq. |6| , we assume that only par- 
ticles bound to the Globular Cluster feel dynamical friction, 
whereas unbound (stripped) particles do not. In this way we 
neglect effects that the density wake self-gravity may induce 
on the trajectories of escaping particles. 

After calculating the different force terms, Superbox 
solves the equation of motion of each cluster particle 

^ = -V$NFw(r0 + ff d£ (n), (7) 

where £ = 1 for bound particles and £ — for unbound ones. 

Superbox uses a leap-frog scheme with a constant 
time-step to integrate the equations of motion for each par- 
ticle. Each N-body realization is evolved a Hubble time 
(t H =14 Gyr) in the NFW dwarf potential. We choose a 
time step dt = t cr /25, where t CI is the cluster crossing time 
defined as 

R 3/2 

The crossing times of Fl, F2 and F3 are, respectively, 15.7, 
3.1 and 0.32 Myr, spanning a factor i C r,Fi/tcr,F3 — 50. Note 
that the crossing time of the Fornax dwarf is t CI ~ 230 Myr 
assuming a dynamical mass of 10 7 Mq and R c = 400 pc, il- 
lustrating the large range of time-scales that our simulations 
must cover, since clusters evolve a factor 10-1000 faster than 
the host galaxy. The total number of time-steps required to 
simulate each of the clusters for a Hubble time is therefore 
large, 1.5 x 10 4 , 8.0 x 10 4 and 7.5 x 10 5 for the clusters Fl, F2 
and F3, respectively. We find that our choice of numerical 
parameters leads to a total energy conservation better than 
5% after the cluster models are evolved for a Hubble time 
in isolation, which yields a negligible evolution of the initial 
model parameters. 



3.2 Relaxation processes in clusters 

Superbox is a collision-less N-body code and therefore does 
not account for a number of internal processes that take 
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place in clusters, such as mass loss due to stellar evolu- 
tion, core collapse and the subsequent re-expansion due to 
two-body relaxation processes. Whereas the evolution asso- 
ciated with stellar evolution is significant only during the 
early life of a cluster ( J; 1 Gyr after formation, Baumgardt 
& Makino 2003), two-(and higher order)-body interactions 
between stars may affect clusters during their entire evolu- 
tion. A realistic incorporation of these effects into dynam- 
ical simulations has proven to be a theoretical as well as a 
numerical challenge due to the complexity of the problem 
and the large range of time-scales involved (e.g. Aarseth & 
Heggie 1998). Only in the last few years has the availabil- 
ity of GRAPE special-purpose supercomputers provided the 
necessary computational tools to simulate the evolution of 
GCs with a number of particles N 10 5 (e.g. Baumgardt & 
Makino 2003). 

To estimate the importance of two-body relaxation pro- 
cesses in the Fornax clusters we calculate first the relaxation 
time as (Binney & Tremaine 2008, eq. 1.38) 



trel 



N-4. 



8 ln(JV*7) 



(9) 



where A r * « Mgc/iti and m is the mean stellar mass in a 
cluster. A Kroupa (2001) mass function then implies m ~ 
0.58M Q . 

Two-body interactions increase the strength of dynam- 
ical friction. As a result, the values of the Coulomb loga- 
rithm in multi-mass clusters are higher than those found 
in collision-less systems. The factor 7 in the Coulomb log- 
arithm depends on the mass spectrum and is rather uncer- 
tain for clusters. As an additional complication, the mass 
spectrum changes during the dynamical evolution of a glob- 
ular cluster (Baumgardt & Makino 2003). Here we adopt 
the value found by Giersz & Heggie (1996), 7 = 0.02. For 
this choice of parameters, we estimate that the Fornax clus- 
ters Fl, F2 and F3 have relaxation times of i re i — 14.6,4.4 
and 2.2 Gyr, respectively. Thus, only Fl has a relaxation 
time trei~te, which suggests that stellar encounters have 
played only a small role over its lifetime. Fornax's other 
clusters, however, have likely suffered repeated episodes of 
core collapse and continuous mass dissolution (Henon 1960). 
In practice, this result implies that their present properties 
have changed significantly since their formation and cannot 
be extrapolated to several Gyr in the past. 



4 MASS AND ORBITAL EVOLUTION OF 
CLUSTERS IN DWARF GALAXIES 

Two parameters govern the dynamical evolution of a gravi- 
tationally bound system within a host galaxy: (i) the mean 
density, which regulates the strength of external tides and 
(ii) the cluster mass, which determines the amount of dy- 
namical friction. 
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Figure 2. Mean density (p*) = Mqc /(<iirR^ /3) versus stellar 
mass for the GC population in the Fornax (filled symbols) and 
Sagittarius (open symbols) dSphs. Massive, low-density clusters, 
which experience the strongest dynamical friction and tidal mass 
stripping, are absent in dwarf galaxies. 



Ter 8 and Fl, will suffer a small or negligible friction force, 
but will respond heavily to tidal interactions with the host 
galaxy. In contrast, owing to its high density and low mass, 
Ter 7 will likely experience little dynamical evolution. At 
high mass and high density, clusters F3, F4, F5 and M54 
may suffer a strong orbital decay as a result of dynamical 
friction. However, since these clusters are considerably more 
dense than the underlying dark matter halo (see FigfTJl, they 
should be resilient to tidal mass stripping. 

Despite the fact that the low-mass end of the For- 
nax/Sagittarius GC distribution shows a remarkable dis- 
persion in stellar densities (~ 2.5 orders of magnitude), no 
massive, low-density clusters are observed in these dSphs. 
The absence of massive, low-density GCs may reflect pro- 
cesses relevant to GC formation or, alternatively, may in- 
dicate that such clusters have been disrupted by tides as 
dynamical friction brought them close to the inner-most re- 
gions of their host galaxies. If the latter scenario is correct, 
the effects of dynamical friction and tidal stripping should 
give rise to spatial segregation according to GC mass and 
density. Specifically we expect surviving low-mass GCs to 
populate the outskirts of their host halos, and we expect the 
most massive GCs to have been dragged to the innermost 
regions, where they can survive only if sufficiently dense. 



4.1 Mass vs density 

In Fig. [2] wc plot cluster mass against the mean density for 
the nine GCs of Fornax and Sagittarius. The large range 
of properties exhibited by these clusters suggests that tidal 
forces and dynamical friction act with different strengths. 
In particular, low-mass, low-density clusters, such as Arp 2, 



4.2 Density vs projected position 

Indeed, Fig. [3] shows that the Globular Clusters of For- 
nax and Sagittarius exhibit a clear spatial segregation. Dots 
show the estimated density of the underlying dark matter 
halos and arrows indicate that the projected radius is the 
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Figure 3. Moan stellar density of clusters (p Jr )/3 as a function 
of their projected position within the dwarf. Arrows indicate that 
the projected position is always smaller than or equal to the dis- 
tance to the dwarf's centre. For comparison, we also plot the 
mean density of the host's halo (pNFw}( r ) f° r the range of max- 
imum velocities estimated for the Milky Way dwarfs. Note that 
the densest and most massive clusters tend to be located in the 
inner regions of dwarfs, suggesting that the cluster population 
may have experienced significant dynamical evolution. Also note 
that the fact that (p*}/3 > (pnfw) suggests that tidal mass loss 
no longer affects these systems at their present location. 
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Figure 4. Mass fraction enclosed within the tidal radius of the 
clusters Fl, F2 and F3 as a function of the distance to the host 
centre. Note that the cluster F3 (also F4 and F5) is so dense that 
it would survive tidal disruption in the potential of any known 
Milky Way dSph. 



minimum possible distance to the dwarf centr^] i.e R < r. 
As expected, Arp 2, Ter 8 and Fl — the lowest-density clus- 
ters in our sample — are found farthest from the dwarf cen- 
tre, whereas M54 and F3 — the densest and most massive 
clusters — reside in the inner regions. 

A second interesting result is that GCs orbiting dwarf 
galaxies tend to reside at a projected radius, 7?, such that 
the mean density of the cluster (p*) = Mgc / (47T-Rf /3) is 
(p*}/3<; (pnfwH-R), which suggests that tidal mass strip- 
ping no longer affects the existing GCs at their present lo- 
cations (see § 4.3 1. As in the case for the mass-density trend 
shown in Fig. [2j it is unclear whether the correlation be- 
tween cluster density and position apparent in Fig. [3] arises 
from a formation mechanism or results from the disruption 
of low-density systems. In this contribution we focus on the 
latter scenario. We begin by examining the effect of dynam- 
ical processes on the evolution of GCs in a dwarf galaxy 
potential. 



4.3 How resilient are GCs to tidal forces? 

We quantify the expected influence of tides by computing 
the theoretical "tidal radius" of the cluster, r t id, as 



(p*>(ftid)/3 = (pnfw)M 



(10) 



where (pnfw) denotes the mean enclosed density within the 
cluster's galactocentric radius r. The estimated amount of 
bound mass is Mgc = Mgc(< rta). In Fig. [4] we plot this 
quantity for the clusters Fl, F2 and F3 as a function of 
their distance to the Fornax centre. Our simple analytic es- 
timates provide a simple means of predicting the dynam- 
ical evolution of GCs. First, the cluster F3 (and similarly 
F4 and F5) is sufficiently dense (see Fig. [T]) that it should 
withstand the tidal field of Fornax even in the inner-most 
regions of the galaxy. In fact, since the Milky Way dwarf 
spheroidals have peak velocities V ma x ~ 40 km/s (Penarrubia 
et al. 2008a,b), our analytical estimates indicate that, if two- 
body relaxation processes are ignored, the clusters F3, F4 
and F5 would likely survive tidal disruption in any of those 
galaxies. Second, intermediate-density clusters, such as F2, 
will suffer tidal mass loss if their orbits bring them to galac- 
tocentric distances similar to their spatial size (r ^ Rt,Gc)- 
Finally, these estimates show that low-density clusters like 
Fl may suffer strong episodes of tidal mass stripping (and 
eventually, disruption) for a large distance range. Stellar de- 
bris resulting from the cluster disruption will dissolve in the 
dwarf galaxy and may be detectable with the present instru- 
mentation. We shall return to this point below. 



4.4 Dynamical evolution of clusters 

To examine the analytical results outlined in §4.3| we have 
carried out a simple numerical experiment. Fig. [5] shows 
the evolution of the galactocentric distance (left panels) and 
bound mass fraction (right panels) for N-body realizations 
of the clusters Fl, F2 and F3 in a dwarf galaxy potential 
that has parameters similar to those of Fornax (see j|2.1[). 



1 M54 sits in the centre of the Sgr dwarf. To plot this cluster in 
Fig. [3] we adopt a galactocentric distance equal to its King tidal 
radius 



© 0000 RAS, MNRAS 000, 000-000 



Tidal disruption of Globular Clusters in Dwarf Galaxies 7 




5 10 5 10 



t (Gyr) t (Gyr) 

Figure 5. Galactocentric distance (left column) and mass (right 
column) evolution of the clusters Fl, F2 and F3 in a dwarf galaxy 
halo with a peak velocity V^nax = 20 km/s. Different colours 
denote different initial galactocentric distances (ro). The orbits 
are confined to the X — Y plane (i.e., the plane formed by the 
halos's long and intermediate axis) and have an orbital circularity 
77 = 1 (i.e. the initial velocity vq is such that fo(ro) = V c (ro)). 



For simplicity, we consider orbits confined to the X-Y plane 
(i.e. the plane formed by the long and intermediate axis of 
the dwarf's halo, respectively). We place model GCs at dif- 
ferent galactocentric distances r(t — 0) = ro with an initial 
velocity that corresponds to the circular velocity of the For- 
nax dwarf at the cluster's initial location, i.e vo = V c (ro). 
The oscillation of the orbital radius around its initial value, 
ro, results from the triaxiality of the halo potential (see e.g. 
§ 3.3 of Binney & Tremaine 2008). 

As expected from the estimates shown in Fig. [4] tidal 
stripping has a minor impact on the evolution of the clus- 
ters F2 and F3. Neglecting the internal evolution induced by 
relaxation processes in these systems, our simulations show 
that the tidal field of Fornax can strip a maximum 20% and 
10% of the initial mass of the clusters F2 and F3, respec- 
tively. 

In contrast dynamical friction induces a strong orbital 
decay of F2 and F3 only for those orbits that are initially 
placed within ro 5; 1.5 kpc from the dwarf galaxy centre. The 
orbital decay process is further strengthened by the small 
amount of mass that these dense clusters lose. The effects 
of dynamical friction become considerably weaker as we in- 
crease ro. F2 and F3 models that are initially placed beyond 
1.5 kpc do not sink to the galaxy centre in a Hubble time. 
These results are in good agreement with those of Goerdt 
ct al. (2006), considering that these authors use a different 
numerical setup. Their models assume point-mass clusters 
orbiting in a live, self-consistent NFW halo with parameters 
similar to those adopted here. The fact that we derive similar 
decay times lends support to the semi-analytical treatment 
of dynamical friction presented in f|3] 

Although the present distances of the clusters F2, F3, 



F4 and F5 to the centre of Fornax are unknown, their pro- 
jected locations fall within the dwarf's limiting radius (see 
Table [TJ . If we adopt the projected position as the galacto- 
centric distance, the fact that these clusters did not sink to 
the galaxy centre within a Hubble time suggests that they 
formed in the outer regions of Fornax and that dynamical 
friction has had little influence on their orbital evolution. 
This explanation, however, places strong constraints on the- 
oretical models of the formation of GCs in dwarf galaxies, a 
topic that still lacks theoretical underpinning. 

We note here that Hernandez & Gilmore (1998) and 
Goerdt et al. (2006) propose an alternative explanation for 
the non-central positions of the Fornax GCs. These authors 
show that adopting a cored halo profile can efficiently stall 
the orbital decay process. Numerical simulations show that 
in cored haloes dynamical friction is suppressed within the 
halo core (Read et al. 2006) and that, as a result, the orbital 
decay of GCs stalls at the core radius and does not proceed 
to the galaxy centre. However, M54 presents a challenge for 
this solution. M54 is the densest and most massive cluster 
found in the Milky Way dSph population and presently sits 
in the centre of the Sagittarius dSph. Its chemical compo- 
sition and kinematics are clearly distinct from those of the 
surrounding stars (Bellazzini et al. 2008), which indicates 
that M54 formed elsewhere and was dragged to its present 
location by dynamical friction, in a fashion similar to our 
models shown here. 

Finally, Fig. [4] shows that the cluster Fl suffers a neg- 
ligible orbital decay regardless of the initial galactocentric 
radius ro. Tidal mass stripping has negligible impact on the 
models that orbit beyond Si 0.5 kpc. Smaller galactocentric 
distances, however, lead to the complete tidal disruption of 
Fl within a Hubble time. In the following Sections, we use 
Fl as a testbed to study the observational imprints that the 
tidal disruption of globular clusters may have left on dSph 
stellar populations. 



5 DISRUPTION OF CLUSTERS IN TRIAXIAL 
POTENTIALS 

5.1 The orbits 

According to CDM modelling, dwarf galaxies are embedded 
in dark matter halos that have a triaxial shape. Orbits in tri- 
axial potentials are usually classified in four major families: 
(i) boxes, (ii) loops (also called tubes), (iii) resonant orbits 
and (iv) irregular (also called unstable or stochastic) or- 
bits. The fraction of phase-space volume occupied by each of 
those families depends on several factors, for example mass 
profile and shape of the galaxy (e.g. Schwarzschild 1979, 
1982; de Zeeuw & Merritt 1983; de Zeeuw 1985; Miralda- 
Escude & Schwarzschild 1989; Merritt & Valluri 1999). The 
study of orbits in triaxial potentials is a complex topic and 
beyond the scope of this contribution. We refer interested 
readers to § 3.3 of Binney & Tremaine (2008) and the pa- 
pers listed above for more detailed analysis. 

For simplicity, we restrict our study to orbits that ini- 
tially lie on the X — Y and X — Z symmetry planes; recall 
that the coordinates (X, Y, Z) are oriented along the ma- 
jor, intermediate and minor axis, respectively. Our N-body 
models are placed along the major axis at an initial distance 
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Figure 6. (X, U) surface of section of the orbits placed at an 
initial distance ro = 0.5 kpc from the dwarf centre. The orbital 
planes are defined by the major and intermediate axis (lower 
panel) and by the major and minor axis (top panel). The pa- 
rameters of our dwarf galaxy halo have been scaled to the values 
estimated for the Fornax dSph, V ma x = 20.6 km/s and r max = 4.1 
kpc. 



ro. The range of initial distances explored here, ro € [0,3] 
kpc, extends considerably farther than the stellar truncation 
radius of Fornax (Rt,For ~ 2.1 kpc). To select the initial ve- 
locity, vo, we define the orbital circularity rj = J/J C (E), 
where J = rovo is the angular momentum and J C (E) the 
angular momentum of a circular orbit with the same energy 
E — 1/2vq + $NFw( r o, 0, 0). This parameter has values in 
the range rj 6 [0, 1]. Because in triaxial systems the angu- 
lar momentum is not a constant of motion, rj is only well 
defined at t = 0. However, introducing this parameter to 
define our initial set of orbits allow us: (i) to vary the initial 
orbital energy at a fixed orbital apocentre in a simple way 
and (ii) to facilitate a comparison with models in spherical 
potentials in order to address the effects of triaxiality (see 



S6.ll 



Fig. [6] shows the surface of section of the orbits placed 
initially at ro = 0.5 kpc. In this plot, each dot corresponds to 
a crossing of the plane specified by Z = Y = 0. The first re- 
sult that stands out is the chaotic nature of the orbits, in the 
sense that small variations in the initial velocity/circularity 
can lead to completely different orbit configurations. This 
plot includes examples of the four orbit families mentioned 
above 

• Box orbits: Systems on box orbits repeatedly pass 
close to the centre of the potential, and eventually at arbi- 
trarily close distances. They have no net angular momentum 
and therefore no sense of rotation. Examples are the orbits 
with rj < 0.2 in both panels of Fig. [6] 

• Resonant orbits: Resonant orbits satisfy a relation of 
the form Iwx +mWY +nwz = 0, where Wi (i = X, Y, Z) are 
the fundamental frequencies around the principal symmetry 




Figure 7. Orbits on the X — Y (upper panel) X — Z (lower panel) 
planes with orbital apocentre ro = 0.5 kpc for different values of 
circularity rj. 



axis and I, m, n integers, not all of which are zero. In the 
surface of section plotted in Fig. [SJ resonant orbits can be 
seen as chains of islands (for example, orbits with rj = 0.3 in 
both panels. Also orbits with rj = 0.5, 0.6 and 0.7 moving on 
the X — Z plane are clear examples of resonant orbits). An 
aspect of these orbits that will prove crucial for our study is 
that, in contrast to boxes, resonant orbits are centrophobic, 
i.e., they avoid the centre of the potential. 

• Loop orbits: Systems moving on loops have a definite 
sense of rotation about the galaxy centre and also avoid the 
region of the origin. Interestingly, we find loop orbits only 
on the plane X — Y. On the plane X — Z the region of 
phase-space occupied by these orbits is filled with resonant 
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and unstable orbits. As resonant orbits, loops are also cen- 
trophobic and do not approach the galaxy centre. 

• Unstable orbits: Unstable orbits have only two con- 
stants of motion in a three-dimensional space. As a result, 
they do not show a definite structure in phase-space (see the 
orbit 77 = f on the upper panel of Fig.j6|. In contrast to box, 
loop and resonant orbits, systems moving on unstable orbits 
will quickly move off the symmetry plane^] These orbits can 
be found more frequently around the minor and intermediate 
axis and do not exist on the plane formed by the major and 
intermediate axis (e.g. Goodman & Schwarzschild 1981). 

We emphasize that structure of the surface of section 
strongly depends on the apocentric radius of the orbit r a « 
ro. As Gerhard & Binney (1985) show, the fraction of loop 
orbits tends to increase at large distances, whereas resonant 
orbits dominate the phase space volume close to the galaxy 
centre. 

The upper and lower panels of Fig. [jj show the projec- 
tion of the orbits examined above on the orbital planes X — Y 
and X — Z, respectively. These Figures illustrate the richness 
of the orbit zoo allowed by a triaxial potential, even if this is 
just a small sample of all the possible orbital configurations. 

A clear transition between box and loop orbits can be 
seen in the upper panel of Fig. [7J as we increase the value of 
the orbital circularity. Orbits with r\ > 0.6 are loops, whereas 
orbits with r\ — 0.1,0.2,0.4 and 0.5 are boxes and boxlets. 
The only example of a resonant orbit on this plane occurs 
for r\ — 0.3, which corresponds to a resonance of the type 
(l,m,n) — (—5,4,0). The number of resonances increases 
on the plane X — Z, as we can see in the lower panel of 
Fig. [7] For example, 77 = 0.3 corresponds to a (—5,0,4) 
resonance whereas 77 = 0.6 is a" fish" orbit with a frequency 
ratio of (3,0,-2). Our potential allows for a much larger 
number of resonances, many of which appear for orbits that 
are not initially confined on the symmetry planes (see Merrit 
& Valluri 1999 for details). 



5.2 Cluster survival time 

For a cluster with a fixed density profile, two factors deter- 
mine the amount of tidal mass loss during a Hubble time: 
(i) the pericentric distance of its orbit and (ii) the num- 
ber of pericentric interactions. The orbital family impacts 
both factors, determining how often and how close a system 
can get to the galaxy centre. A cluster moving on a box 
orbit passes repeatedly through the centre of the galaxy, 
whereas a loop orbit always avoids the inner-most regions of 
the potential. As Fig. [6] shows, resonant orbits represent an 
intermediate case. We therefore expect that clusters on box 
orbits will suffer the strongest episodes of mass loss while 
those on loop orbits will have the largest survival times. On 
the other hand, the number of pericentre passages within 
a Hubble time depends on the orbital apocentre. Clusters 
with a large ro spend a large fraction of the orbital time 
at large distances from the galaxy centre and thus suffer a 



2 For a recent study of orbit instability in a triaxial potential such 
as the one considered here see Adams et al. (2007) and references 
therein. 
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Figure 8. Survival time (t s ) of the cluster Fl as a function of the 
initial distance ro to the galaxy centre and the orbital circularity 
rj. Horizontal dotted lines indicate the value of a Hubble time 
tff = 14 Gyr. The NFW potential has a peak velocity V max = 20 
km/s, similar to the value estimated for the Fornax dSph. 



small number of pericentric interactions. As a result, we ex- 
pect that the mass loss rate will decrease on average as ro 
increases. 

We define the cluster survival time, t s , as the time re- 
quired for a cluster to lose 95% of its initial mass. Although 
arbitrary, this definition provides an accurate estimate of 
the survival time of cored clusters. Penarrubia et al. (2002) 
and Zhao (2004) find that cored stellar systems that are 
subject to tidal mass stripping show a characteristic mass 
evolution that can be divided into two different regimes: (i) 
a prolonged and slow mass loss regime that is typically fol- 
lowed by (ii) a sharp tidal disruption that merely spans a 
few crossing times. This implies that our particular defini- 
tion of "survival time" has an inaccuracy of the order of the 
cluster's dynamical time, i.e. a small fraction of a Hubble 
time for the cluster models considered here. 

We must note that, since relaxation processes speed up 
the rate of tidal mass stripping (see j 



3.2 for details), the sur- 



vival times derived here represent upper values. On the other 
hand, the fact that our potential is kept fixed during the 
whole cluster evolution tends to underestimate the survival 
time of GCs in dwarf galaxies because, in self-consisent sys- 
tems, the inner regions of dark matter haloes are expected 
to react when the enclosed halo mass compares with the 
mass of the cluster, which in practice decreases the strength 
of the dwarf's tidal field. 

Bearing these limitations in mind, Fig.[8]shows the sur- 
vival times computed for the cluster Fl in a potential scaled 
to the parameters of the Fornax dSph. We have explored a 
wide range of apocentric distances < ro < 2.5 kpc, and 
note that the stellar truncation radius of Fornax occurs at 
Rt.For = 2.1 kpc. As the Figure shows, the survival time in- 
creases monotonically as the apocentric distance increases. 
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Survival time also tends to increase with orbital circu- 
larity, although the variation is not entirely monotonia The 
oscillations of t s = t s (rf) for a fixed value of ro are visible 
in both panels and arise from changes of the orbital family. 
For example, inspection of the systems with ro = 2.0 kpc in 
the upper panel reveals that (i) values r/ = 0.1,0.4,0.7 cor- 
respond to box orbits and lead to local minima of t s , as one 
would expect, (ii) values r\ — 0.2, 0.3, 0.5 and 0.6 correspond 
to resonant orbits, and (iii) r\ > 0.8 are irregular box orbits. 
The results are similar for other values of ro, although ir- 
regular orbits with high circularity become more common 
as the orbital apocentre increases. 

The strong sensitivity of t s to the orbital family is also 
evident in the lower panel of this Figure. For example, the 
sharp increase in t s for r\ > 0.6, independent of the value 
of ro, corresponds to the transition between box and loop 
orbits (see Fig. [6|. The latter maximize the survival times 
of the cluster and Fig. [8] shows that Fl does not disrupt if 
the loop orbit has r a ^ 1.5 kpc. We expect therefore that the 
cluster debris associated with this orbital family is deposited 
predominately in the inner regions of dwarf galaxies (see 
§[63l. 

We note that a large fraction of box and resonant orbits 
(rj < 0.4) on both orbital planes leads to the tidal disruption 
of the cluster Fl within a Hubble time for the whole range 
of orbital apocentres explored here. Some of these orbits 
bring globular clusters beyond the stellar truncation radius 
of Fornax which, as we will see below, may (possitively) 
impact our ability to detect the debris. 

5.3 Debris morphology 

Stars that escape from a system undergoing tidal disrup- 
tion have energies and angular momenta that are similar to 
those of the progenitor system (e.g. Penarrubia et al. 2006). 
The scatter in these quantities causes a progressive phase- 
mix of debris along the progenitor's orbit until, eventually, 
the full phase-space volume is filled (Johnston et al. 1998, 
Helmi & White 1999). The time that tidally stripped stars 
need to fill the phase space volume of the progenitor's or- 
bit is of the order of the dwarf's crossing time. The fact 
that typical dwarf galaxies have crossing times that are con- 
siderably shorter than a Hubble time (for example in the 
case of Fornax, which has a crossing time t CI ~ 230 Myr, a 
GC may perform of the order of tu /t C i — 60 revolutions) 
suggests that the debris from tidally disrupted GCs fully 
phase-space mix on a relatively short time scale. Also, any 
transient feature is necessarily short-lived and disappears 
within few crossing times. 

Because regular orbits conserve three integrals of mo- 
tion in a static three-dimensional potential, they distribute 
in toroidal surfaces in a six-dimensional (6D) phase space 
(e.g. Binney & Tremaine 2008). Unfortunately, 6D orbital 
information is barren to us, as only projected quantities can 
be known with the present instrumentation. As the catastro- 
phe theory explains (e.g. Arnold 1992), the projection of a 
six-dimensional manifold onto a two-dimensional space nec- 
essarily results in the formation of folds and cusps (see e.g. 
Cappellari et al. 2004 for a schematic representation). Folds 
are curves on the plane that arise as two inverse images of 
an orbit merge and disappear, whereas cusps appear at the 
intersection of two or more folds. Both folds and cusps will 
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Figure 11. Mean line-of-sight velocity of the cluster Fl debris. 
Only models moving on loop orbits on the X-Y plane are con- 
sidered. The line-of-sight projection is aligned with the Y-axis 
to maximize the rotational velocity component. Note that, al- 
though the rotation signal strengthens as the orbital pericentre 
increases, the number of loop orbits that lead to the disruption 
of Fl decreases. We find that maximum rotational velocity can 
reach values of the order of the dwarf's peak velocity (V ma x = 20.6 
km/s). At large radii, however, the mean rotational velocity falls 
to (vi. . s ) ~ 5 km/s. 

appear as prominent features in the projection of cluster 
debris on the sky. 

In Fig. [9] and [To] we show the configurations of tidal 
debris that result from fully disrupted Fl models in a triax- 
ial potential. The orbits of the progenitor systems are those 
shown in Fig. [7] Note that this selection only corresponds to 
a small subset of the whole orbital family spectrum allowed 
by the potential of eq. |4j and that other resonant families 
yield fairly different debris morphologies. The orientation of 
the line-of-sight direction with respect to the progenitor's 
orbital plane also alters the debris morphology. For simplic- 
ity, we only show here the projection along the three main 
orbital axes. Values of the projected surface brightness are 
colour-coded for clarity, and in this Figure bright areas cor- 
respond to regions densely populated by debris material in 
a particular projection. 

These Figures show that the cluster debris neatly traces 
the orbits plotted in Fig. [7] and that, as a result, the pro- 
jected morphologies are extremely rich. Folds and cusps ap- 
pear in practically all orbital configuration and may produce 
striking features. For example, when looking from a direction 
perpendicular to the orbital plane, folds associated with box 
orbits exhibit a "butterfly" -like shape, whereas loops tend to 
generate disc-like configurations of debris with a noticeable 
scarcity of stars in the inner-most regions. Complex stellar 
substructures may arise from the disruption of clusters on 
resonant orbits. A good example is the "fish" orbit shown 
by Fig. |10| Note, however, the strong dependence of debris 
morphology on the line-of-sight projection: (i) its complex- 



© 0000 RAS, MNRAS 000, 000-000 



Tidal disruption of Globular Clusters in Dwarf Galaxies 11 




Figure 9. Projected surface brightness associated with the tidal debris from the cluster Fl on different orbits initially confined on the 
symmetry plane X — Y. The initial apocentric distance is ro = 0.5 kpc and the pixel resolution is 0.04ro- Note that the distribution of 
debris reflects the progenitor's orbit shown in the upper panel of Fig. [7] 



ity decreases if the projection angle is not aligned with the 
axis perpendicular to the orbital plane and (ii) the surface 
brightness and the area occupied by tidal debris decrease if 
the projection angle is aligned with the main orbital axis, 
which in this particular coordinate system corresponds to 
the X axis. 

Our models show that stellar substructures arise nat- 
urally from the disruption of GCs over a Hubble time in 
a static potential. Shells (folds) and isolated over-densities 
(cusps) arise in almost all orbital configurations and line- 
of-sight projections. These features result from the triaxial 
shape of the galaxy potential and would not form in spher- 
ical haloes. Their surface brightnesses range from 36 to 26 
mag/arcsec 2 for a progenitor cluster with an initial lumi- 
nosity of L ~ 10 4 Lq. Note, however, that also clusters with 
higher luminosities can be fully disrupted in the tidal field 
of a dwarf if their initial mean density is (p*) 0.1Z/©/pc 3 
(see f4.3l and that these systems may leave behind stel- 



hinder or even completely prevent their detection in pho- 
tometric surveys. We examine this issue in more detail in 



lar debris with surface brightneses below 26 mag/arcsec 2 . 
Although the latter value is around six magnitudes brighter 
than the typical detection threshold that present instrumen- 
tation can reach for a resolved stellar population, the fact 
that they might be surrounded by a sea of field stars may 



5.4 Kinematics 

The disruption of a GC in a dwarf galaxy also leaves kine- 
matic imprints. The kinematic properties of the debris are 
fairly sensitive to the orbit of the progenitor cluster. For 
example, box and resonant orbits have zero net angular mo- 
mentum, causing the cluster debris to have an average line- 
of-sight velocity that is null everywhere regardless of viewing 
angle. However, loop orbits have a net rotational velocity 
component that will be reflected by the kinematics of the 
tidal debris unless the viewing angle is perpendicular to the 
orbital plane. 

For example, in Fig.[TT]we plot the variation of the mean 
line-of-sight velocity (vi, , s ) for a projection angle parallel to 
the orbital plane. This Figure includes all our Fl models on 
loop orbits that tidally disrupt within a Hubble time. Note 
that no loop orbits lead to the full disruption of the cluster 
Fl for orbital apocentres ro <; 1.5 kpc (see models with cir- 
cularity r) > 0.6 in the lower panel of Fig. RjJ. Models with 
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Figure 10. As Fig.|9]for orbits initially on the X — Z plane 



large apocentres tend to exhibit a high line-of-sight peak ve- 
locity that can reach values comparable to the halo's peak 
velocity, V m ax,Foi = 20.6 km/s. In all the models examined 
here the maximum of {vi. .„) occurs at ~ 0.3-0.5 kpc from 
the dwarf's centre. Beyond that radius this quantity rapidly 
decreases, leveling off at ~ 5 km/s. 

Although loop orbits are interesting due to the rota- 
tional component that they introduce in non-rotating sys- 
tems like dSphs, the dominant orbital families in the inner 
regions of triaxial NFW dark matter halos are box and res- 
onant orbits (see Fig. |6j. To examine the debris kinematics 
associated with these types of orbits we inspect orbits on 
the X-Z plane. 

In Fig. [12] we plot the projected velocity dispersion de- 
fined as uf o s = {(»(.„., — (ui.o.s)) 2 } for the models shown in 
Fig. |10| For comparison purposes, we note that the field stars 
of Fornax show a flat velocity dispersion profile ap r — 10.5 
km/s (Walker et al. 2006a, 2007). 

Fig. [I2]displays a strong dependence of velocity disper- 
sion on the projection angle: cluster debris is fairly cold, 
o";.o.s~3 km/s for all orbital circularities, so long as the 
projection angle is perpendicular to the orbital plane. This 
reflects the low velocity dispersion of the progenitor clus- 
ter Fl <7 = [4ttG/9R 2 c po] 1/2 ~ 1.34 km/s (see Table [lj. In 



this particular line-of-sight projection, all the bright features 
visible in the X-Z panels of Fig. [10] and [12] have a velocity 
dispersion comparable to that of the progenitor cluster and 
always exhibit a high degree of symmetry with respect to the 
dwarf centre regardless of the progenitor's orbital family. 



Contrary to naive expectations, owing to the high radial 
anisotropy of these orbits the cluster debris can also appear 
hotter than the surrounding field stars if the line-of-sight 
direction lies on the orbital plane. As Fig. [12] illustrates, 
cluster debris may have projected velocity dispersions with 
strong gradients throughout the galaxy and a fairly complex 
geometry. Box orbits tend to show a decreasing dispersion 
profile from the dwarf's centre, whereas off-centre, isolated 
hot clumps arise from resonant orbits. In these models, there 
is a clear correspondence between the brightest areas shown 
in Fig. [TO] and the "hottest" regions visible in Fig. |12| which 
show projected velocity dispersions that can be as high as 
~ 14 km/s. 
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Figure 12. Projected velocity dispersion of the models shown in Fig. |10| Note that all models shwon here move on either box or resonant 
orbits. For comparison, Fornax's field stars have a projected velocity dispersion of ~ 10.5 km/s practically constant at all radii (Walker 
et al. 2007). Note the strong dependence on the line-of-sight projection: debris appear fairly cold (c;. . s < 3 km/s) if the line of sight 
is perpendicular to the orbital plane. However, if the line of sight is aligned with the orbital plane, strong velocity dispersion gradients 
arise in almost all orbital configurations. In localized areas the tidal debris field may appear considerably /io44er than Fornax's field stars, 
with a projected velocity dispersion that may reach <X(. o s ~ 14 km/s. 



6 DISCUSSION 

6.1 Effects of halo triaxiality 

As previously discussed, a fundamental prediction of the 
CDM paradigm is that dark matter haloes have a triax- 
ial shape. There are few contributions in the literature that 
study the disruption of stellar systems in triaxial potentials, 
and it is interesting to analyze here, though briefly, how the 
adoption of a triaxial halo may have influenced our results. 

In contrast to spherical haloes, which allow only rosette 
orbits, triaxial haloes allow a wide range of orbital families 
with very different properties, as outlined in §5.1| In the 
context of the survival of stellar clusters in dwarf galaxies, 
this has strong consequences. 

As loop orbits, rosettes are centrophobic and avoid the 
centre of the potential, with the single exception of a ra- 
dial (r) = 0) orbit. Further, the pericentric distance and the 
orbital period are increasing functions of 77, which implies 
that, at a fixed apocentre, the survival time of star clusters 
t s — t B (rf) increases monotonically. As we discuss in { 5.2 



this characteristic is not shared by box and resonant orbits, 
which can come arbitrarily close to the centre of the host 
in repeated occasions with independence of the value of rj. 
As a result, clusters that move on box and resonant orbits 
have shorter survival times than those on rosette orbits for 
a given orbital circularity and apocentric distance ro. 

To illustrate quantitatively the change in t s introduced 



by the halo triaxiality, we plot in Fig. 13 the ratio t^/t 1 / 1 
of the cluster Fl as a function of circularity for different 
orbital apocentres. This Figure shows that, independently 
of the orbital parameters, triaxiality always strengthens the 
tidal disruption of stellar clusters, so that t® p The 
difference is smallest for the lowest values of 77 (i.e. highly 
radial orbits in spherical potentials) and, due to the simili- 
tude between loop and rosette orbits, for the combinations 
of ro and 77 that correspond to loop orbits in a triaxial po- 
tential (e.g. ro = 0.5 kpc and rj > 0.6 in the lower panel). As 
expected, clusters on box and resonant orbits survive shorter 
times than on rosettes. This result is particularly evident in 
the upper panel of Fig. |13| recalling that in our triaxial halo 
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Figure 13. Ratio between the survival time of the cluster Fl 
in a spherical (t| ph ) and a triaxial (t 1 / 1 ) halo with axis-ratios 
(a, b, c) = (1, 47, 1.22, 0.74). This quantity is plotted as a function 
of the initial distance ro to the galaxy centre and the orbital 
circularity r\. The NFW potential has a peak velocity V ma x = 20 
km/s, similar to the value estimated for the Fornax dSph. 



model all orbits moving on the X — Z plane are either box 
or resonant orbits. The difference in survival time becomes 
stronger as we move to high 77, which is equivalent to in- 
creasing the orbital pericentre in a spherical potential, and 
can be as large as t s s ph ~ 10ts rl for 77 ~ 1 (i.e. circular orbits 
in a spherical potential). 

Loop orbits and rosettes share a good number of char- 
acteristics. For example, not only do both avoid the centre 
of the potential, also both have a definite sense of motion 
around the galaxy centre. The spatial distribution of tidal 
debris from disrupted clusters moving on rosettes will be 
therefore similar to that shown in Fig. [9] for the models with 
77 > 0.6. These tidal debris will also have a mean rotational 
velocity that will extend approximately out to the orbital 
apocentre of the progenitor cluster, in a similar fashion as 
discussed in §5.4| Finally, a noteworthy corollary is that the 
stellar features associated with the tidal disruption of a stel- 
lar cluster moving on box and resonant orbits (e.g. shells and 
isolated clumps, see Fig. [9] and 10 1 do not arise in a spheri- 
cal potential and may be used therefore to address the shape 
of dark matter haloes. 



6.2 Detectability of Globular Cluster debris 



In [ 5.3 we show that the tidal disruption of a W 4 Lq globular 
cluster in the triaxial potential of a typical dwarf galaxy may 
form stellar substructures that do not dissolve in time. With 
surface brightnesses fi <; 26 mag/arcsec 2 , some of these fea- 
tures fall within the detection limits of present instrumen- 
tation. However, most of those substructures will be sur- 
rounded by a sea of dwarf field stars, hindering the iden- 
tification of cluster debris as stellar over-densities in both 



photometric and kinematical surveys. The detection of a dis- 
ruption event will therefore be more likely in regions where 
the surface brightness of the cluster debris is similar to or 
higher than that of the field stars. 

To explore this issue in some detail, Fig. [14] shows the 
surface density profiles averaged on circular annuli of the 
cluster Fl models that result in a complete tidal disrup- 
tion within a Hubble time. We distinguish here between or- 
bits initially confined to the X-Y (left panel) and X-Z (right 
panel) planes. For simplicity we consider a line-of-sight pro- 
jection perpendicular to the orbital plane. We also plot the 
King model fits of the surface density profiles of the Fornax 
and Carina dwarfs (Mateo 1998) with dotted and dashed 
lines, respectively. 

We begin by discussing the orbits with the smallest or- 
bital apocentre considered in this paper, ro = 0.5 kpc. As 
we show in Fig. [8] these orbits result in the complete tidal 
disruption of Fl regardless of its orbital circularity 77. As the 
top-right and left panels show, the surface density profile of 
the tidal debris is sensitive to the orbital family: (i) box (e.g. 
77 = 0.1) and resonant (e.g. 77 = 0.3) orbits yield debris distri- 
butions resembling King (1962) models, i.e. cored in the in- 
ner regions with a central surface density Eo ^ 1O 5 M0 /kpc 2 , 
and having a sharp truncation radius at large radii, (ii) Loop 
orbits (77 > 0.6 in the left panel), in contrast, yield debris 
distributions with a clear scarcity of stars in the inner-most 
regions of the dwarf. In all orbital configurations cluster de- 
bris shows a truncation radius that closely correlates with 
the progenitor's apocentre. In particular, we find that the 
distribution of debris is truncated at R w 1.2ro regardless 
of the orbital family. 

As discussed in §5.2[ as we increase the orbital apocen- 
tre of the cluster Fl, only box and resonant orbits lead to 
a complete tidal disruption within a Hubble time. This is 
because these orbits bring the cluster Fl repeatedly to the 
centre of the dwarf galaxy regardless of the orbital apocen- 
tre. A crucial result for the possible detectability of the tidal 
debris is that disruption occurs even for orbits with apoc- 
entres that are larger than the limiting radius of the host 
galaxy, i.e. ro Si i?t,F or = 2.1 kpc for the Fornax dwarf. As 
the orbital apocentre increases, the cluster debris extends 
over larger volumes, leading to a drop in central surface 
density. We find that for apocentres ro = 1.5 and 2.5 kpc, 
the central surface density of the debris lies in the range 
Eo ~ lO 3 -lO 4 M0/kpc 2 depending on circularity and orbital 
orientation. 

The maximum surface density (averaged over circular 
annuli) associated with tidal debris from the cluster Fl is 
E ~ 10 5 Mq /kpc 2 . This value is considerably lower than that 
of both Fornax and Carina, which implies that identifying 
cluster stars from the background stellar population in the 
central regions of dwarfs will require a sizable data set and 
additional kinematical information. 

In contrast, cluster debris may dominate the overall 
light distribution at large radii if the progenitor cluster was 
formed beyond the limiting radius of the host galaxy (see 
JO) . Unfortunately, Fig. 



14 



shows that even in regions with 
a scarce presence of dwarf stars the detection of cluster de- 
bris may be observationally challenging in extended, lumi- 
nous galaxies like Fornax given the depth that the survey 
should reach, approximately four decades below the central 
surface density. 
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Figure 14. Surface density profiles averaged on circular annuli of the cluster Fl debris. The left and right panels include models that 
lead to the complete tidal disruption of the cluster Fl moving on the X-Y and X-Z planes, respectively. Rows show orbits with different 
orbital apocentres (ro). Dotted and dashed lines show the surface density profiles of Fornax and Carina, respectively. Note that in some 
of these models cluster debris dominates the overall light profile beyond the limiting radius of the dwarf. 



Finally, a noteworthy remark is that sizable spectro- 
scopic surveys may reveal the presence of cluster debris re- 
gardless of the host background surface density if the now- 
defunct clusters had a chemical abundance pattern different 
from that of most of the dSph stars. We shall return to this 
issue below. 

6.3 Were Globular Clusters more common in 
dSphs? 

At present, only the Sagittarius and Fornax dwarfs have 
associated populations of bound GCs. An understanding of 
the formation mechanism of these objects is still lacking, 
and it is unclear whether other dwarf galaxies may have also 
formed clusters that were subsequently tidally disrupted. It 
is thus interesting to carry out a Gedankenexpervment that 
explores whether the set of observables in the present dSph 



population could help us to infer if such disruption events 
occurred in the past. 

As discussed in the previous Section, the fact that For- 
nax and Sagittarius are the most luminous and spatially 
extended dSphs of the Milky Way (Mateo 1998) may hide 
possible underlying stellar over-densities with low surface 
brightness. Fig. |14| shows that the potential presence of GC 
debris would be considerably more straightforward to de- 
tect in small, faint dwarf galaxies like Carina. According to 
the estimates of Penarrubia et al. (2008a), Fornax and Ca- 
rina are embedded in dark matter halos with similar NFW 
parameters, despite the fact that Fornax is a factor ~ 100 
more luminous and has a half-light radius that is around 
three times larger. The similarity of their halos implies that, 
for a given orbit, the cluster Fl would follow practically the 
same dynamical evolution in both galaxies. Taking the Ca- 
rina dwarf as an example, the tidal disruption of the cluster 
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Figure 15. Evolution of the spatial and kinematical distribution of cluster Fl debris in an NFW halo that loses a fraction of its initial 
mass to tides. This particular model corresponds to the orbit with ro = 0.5 kpc, r\ = 0.8 shown in Fig. |10| and |12| Note that the tidal 
stripping of dwarf galaxies may potentially erase stellar substructures resulting from the tidal disruption of stellar clusters. 



Fl would be easily detectable at approximately two decades 
below the central surface brightness as an excess of light be- 
yond the King limiting radius derived from the inner profile. 
Interestingly, such excesses have been found in several dwarf 
galaxies (Carina, Majewski et al. 2005; Leo I, Sohn et al. 
2007; Sculptor, Westfall et al. 2006; Ursa Minor, Martmez- 
Delgado et al. 2001, Palma et al. 2003), although they have 
been often interpreted as signatures of a recent tidal inter- 
action with the Milky Way (e.g. Mufioz et al. 2008 and ref- 
erences therein). However, Penarrubia et al. (2008c) argue 
that the orbital motion of some of these systems (e.g. Leo 
I) is incompatible with the presence of unbound stellar ma- 
terial within the area surveyed by current data, which casts 
doubts on the tidal origin of the reported extra-tidal features. 

The models presented in the previous Sections show 
that kinematic information may help to identify globular 
cluster debris. As discussed in §5.3| and |5.4| if the cluster's 
orbital plane is approximately perpendicular to the line-of- 
sight directions, shells and off-centre clumps with cold pro- 
jected velocity dispersions would be a clear signature of a 



cluster merger. Encouragingly, the presence of stellar sub- 
structures has been reported in Fornax (Coleman et al. 2004, 
2005; Olszewski et al. 2006), Canes Venatici (Ibata et al. 
2006, Martin et al. 2008; also see Simon & Geha 2007 for a 
discrepant result) and Sextans (Kleyna et al. 2004, Walker 
et al. 2006b), although some of these detections remain con- 
troversial due to their marginal statistical significance. 

The presence of an off-centre localized stellar clump, 
reminiscent of those detected in dwarf elliptical galaxies (e.g. 
Binggeli et al. 2000), is probably clearest in the Ursa Minor 
dwarf (Olszewski & Aaronson 1985; Irwin & Hatzidmitriou 
1995; Wyse et al. 2002; Palma et al. 2003), and has also been 
confirmed by deep photometric data taken with the Wield 
Field Planetary Camera (WFPC2) on the HST (Battinelli & 
Demers 1999). This over-density locates 13" (~ 0.021 core 
radii) from the dwarf's centre. The metalicity of the stel- 
lar clump appears indistinguishible from the bulk popula- 
tion of the dwarf (Eskridge & Schweitzer 2001), which only 
contains a single, old stellar population (e.g., Mateo 1998; 
Feltzing, Gilmore & Wyse 1999). Interestingly, Kleyna et al. 
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(2003) also reported the presence of a kinematic substruc- 
ture in UMi that locates at a projected distance of 0.23 kpc 
from the dwarf's centre (i.e. 1.2 core radii) and has a very 
low velocity dispersion, ~ 0.5 km/s. The UMi substructures 
detected in photometric and spectroscopic data may corre- 
spond to the debris of one or more disrupted stellar clusters. 
However, it must be noted that in our triaxial halo models 
cold features always show a symmetric pattern with respect 
to the dwarf centre. The fact that only a single kinemati- 
cally cold substructure has been detected may imply that 
a significant fraction of tidal debris remain undetected or 
that the origin of the cold clump is not related to a recent 
merger event whatsoever. A third possibility was proposed 
by Keyna et al. (2003), who showed that dissolved stellar 
clusters in spherical haloes may retain a coherent and cold 
spatial distribution for several orbital periods if the halo pro- 
file is cored. The large photometric and spectroscopic studies 
currently underway will be able to quantify the presence or 
absence of symmetrically-distributed substructures in UMi, 
and provide a clean distinction between the two possible 
histories: tidal disruption in a triaxial cusped potential, and 
quiet evaporation in a cored spherical potential. 

Projection angles that are strongly misaligned with the 
normal vector of the orbital plane are more uncongenial 
to the identification of past mergers and lead to a large 
variety of possible spatial and kinematical configurations 
of debris that go beyond the expected "cold stellar over- 
densities". For example, if the line-of-sight direction lies 
close the progenitor's orbital plane, the associated debris 



may show strong velocity dispersion gradients (see Fig. 121. 
Also, the tidal debris from a globular cluster that moves on 
a loop orbit may introduce a rotational velocity component 
in the stellar population of the dwarf. As discussed above, 
a rotation signal may be straightforward to detect in dwarf 
galaxies, which are prototypical examples of non-rotating, 
pressure-supported systems, especially at radii where the 
overall surface brightness is dominated by cluster stars, e.g. 
RftO.G kpc in the Carina dwarf. Signatures, or hints, of ro- 
tation have been reported in kinematic surveys of Carina 
(Munoz et al. 2006; Walker e tal. 2008), Leo I (Mateo et al. 
2008), Sculptor (Battaglia et al. 2008; Walker et al. 2008) 
and Fornax (Walker et al. 2008). A noteworthy remark is 
that tidally unbound stars escaping the system through tidal 
tails would also produce velocity gradients along dSphs sim- 
ilar to those plotted in Fig. [II] Also the transverse bulk mo- 
tion of the dSph galaxy can cause small velocity gradients 
(Kaplinghat & Strigari 2008; Walker et al. 2008). Further 
data that focus on the rotating stellar component is needed 
to distinguish between these scenarios. 

In this context, chemical tagging can possibly provide 
the most reliable means of shedding light on the nature of 
the stellar over-densities and kinematic oddities detected in 
several dSphs. Recent measurements of the chemical abun- 
dances of different elements of the Fornax globular clusters 
(Letarte et al. 2006, Mottini & Wallerstein 2008) reveal 
a composition that is similar to Galactic globular clusters 
and more metal-poor than that of the Fornax field stars. If 
the same result can be extrapolated to the globular clus- 
ters that may have inhabited other dwarf galaxies, a chemo- 
dynamical mapping of dwarf galaxies would possibly provide 
the most efficient method to address the origin of light ex- 
cesses at large radii. 



6.4 Evolving dwarf galaxy potentials 

The results obtained in the previous Sections assume a static 
dSph potential throughout the evolution of our GC models. 
Although this approximation may be appropriate for iso- 
lated dwarf galaxies, it is unlikely to hold in detail for "sub- 
haloes" orbiting in the main halo of the Milky Way. Recently, 
Stoehr et al. (2002), Hayashi et al. (2003), Kazantzidis et al. 
(2004a) and Diemand, Kuhlen & Madau (2007) have stud- 
ied the modifications suffered by NFW haloes undergoing 
tidal disruption within a more massive system. A recurrent 
result is that stripping occurs predominantly in the outer 
regions of the halo and reduces both Vmux and r max . The 
reduction in r max is larger, so that the halo concentration 
increases. As Hayashi et al. (2003) and Penarrubia et al. 
(2008b) show, the variation of both parameters is solely con- 
trolled by the amount of dark matter lost to tides. Because 
the stellar components of most Milky Way dSphs are deeply 
segregated within the dark matter haloes (Penarrubia et al. 
2008a, b) , dSphs must lose at least 90% of their initial mass 
before starting to shed their stars. 

To illustrate how the stripping of the outer regions of 
dwarf galaxy haloes may affect existing stellar substructures, 
we have re-integrated the models shown in previous Sections 
for an additional Hubble time in dark matter haloes that 
have lost respectively 90% and 99% of their initial halo mass. 
Applying the results of Penarrubia et al. (2008b), the halo 
peak velocity Kiax drops to ~ 0.67 and 0.33 of its initial 
value and the peak velocity radius r max becomes a factor 
0.37 and 0.14 smaller, respectively. 

Figure fl5] shows how the tidal debris of the cluster Fl 
respond to the potential variation. For illustrative purposes 
we choose an orbit that leads to the formation of a substan- 
tial amount of substructure (ro = 0.5 kpc, r\ = 0.8, orbital 
plane X-Z) and assume that tidal mass loss does not alter 
the shape of the halo. 

This Figure illustrates a few interesting points. One is 
that the volume occupied by the cluster tidal debris expands 
as a result of lowering the potential depth. This in turn 
causes an overall decrease of the surface brightness and the 
prominence of the brightest over-densities wanes as a result. 
For example, the two bright clumps visible in the Y — Z 
(upper-right) panel of the unperturbed model (M/Mo = 1) 
disappear when the halo loses 90% of its initial mass (i.e 
M/Mq = 0.1). A similar effect is visible in all line-of-sight 
directions, which clearly indicates that the amount of stel- 
lar substructure decreases if a dwarf galaxy undergoes tidal 
mass stripping. In fact, stronger mass loss episodes may 
erase these stellar over-densities altogether, as the right pan- 
els (M/Mq = 0.01) show. Cluster debris integrated in a halo 
that has lost 99% of its mass appear basically as a smooth 
ellipsoid with no obvious substructure visible in any projec- 
tion angle, despite that the halo potential is indeed triaxial. 
Finally, another interesting point may be gleaned from the 
lower panels by noting that the tidal debris become over- 
all "colder". The strong velocity dispersion gradients that 
are visible in the unperturbed model (left panels) tend to 
smooth out as kinematically distinct clumps reduce their 
size (see bottom-middle panel). As in the surface brightness 
maps, strong episodes of mass loss (i.e M/Mq < 0.01) erase 
any previously existing kinematic substructure in the dwarf. 
These experiments therefore suggest that dSphs tend to be- 
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come spatially and kinematically uniform as they lose their 
outer dark matter halo envelopes to tides. 



7 SUMMARY 

We have used N-body simulations to study the tidal dis- 
ruption of Globular Clusters (GCs) in dwarf spheroidals 
(dSphs). We assume a cosmologically motivated scenario 
wherein dSphs are dark-matter dominated systems that fol- 
low an NFW profile and have a triaxial shape. GCs are mod- 
elled as N-body realizations of a King (1962) sphere. We use 
the Fornax and Sagittarius dwarfs, which contain five and 
four GCs respectively, as testbeds for analyzing the orbital 
and mass evolution of their GC populations under different 
orbital configurations. Our findings may be summarized as 
follows: 

(i) GCs in dSphs span over three orders of magnitude in 
density and two orders of magnitude in mass, in contrast 
with the narrow density range estimated for the inner re- 
gions of the dark matter halos in which dwarf spheroidals 
are embedded. For example, the value of the halo density 
measured at 100 pc from the halo centre only varies a factor 
~ 3 among all the Milky Way dwarfs. This result implies 
that the dynamical evolution of the Fornax clusters can be 
safely generalised to the rest of dwarf galaxies that contain, 
or used to contain, Globular Clusters. 

(ii) Clusters in dSphs show a remarkable spatial segre- 
gation depending on their mass and density, such that the 
densest and most massive clusters are located at close (pro- 
jected) distances from the dwarf. This may result from or- 
bital decay induced by dynamical friction. 

(iii) We show that the most massive clusters in For- 
nax and Sagittarius (Mac ft 1O 5 M0) have central densities 
(p*(0) >10 3 M Q /pc 3 ) that are orders of magnitude higher 
than those of the underlying dark matter halos. Our N-body 
simulations show that these systems will survive tidal dis- 
ruption in any dwarf galaxy of the Milky Way. 

(iv) In contrast, clusters with masses < 10 4 M Q , like the 
Fornax cluster Fl, are much less dense (p*(0) 3Mq/pc 3 ) 
and can be fully tidally disrupted if the orbits bring them 
close to the dwarf centre. That introduces a large depen- 
dence between the cluster's survival time and the orbital 
type. At a fixed orbital apocentre, we find that box, reso- 
nant and loop orbits lead to the shortest, intermediate and 
longest survival times, respectively. 

(v) The tidal disruption of a 1O 4 L0 cluster results in the 
formation of stellar substructures that have surfaces bright- 
nesses fi <^ 26 mag/ arcsec 2 and very complex spatial mor- 
phologies. Depending on the cluster's orbital family and line- 
of-sight projection, we find that cluster debris may appear 
as shells, isolated clumps and elongated over-densities. 

(vi) The kinematics of the cluster debris are strongly sen- 
sitive to the relative orientation between the line-of-sight 
direction and the cluster's orbital plane. On average, per- 
pendicular and parallel orientations yield substructures that 
have, respectively, lower and higher projected velocity dis- 
persions than the field stars of the host dwarf galaxy. 

(vii) We find that the disruption of GCs moving on loop 
orbits may introduce a rotational component in dSphs. Our 
models show that, depending on the orbital apocentre, tidal 
debris may reach a maximum mean rotational velocity that 



is comparable to the peak velocity (I^nax) of the host galaxy 
halo. 

(viii) A crucial remark is that these substructures arise 
from the potential triaxiality and do not have a transient 
nature. Indeed, the presence of distinct stellar substructures 
in dSphs may be a signature that gravitationally bound sys- 
tems were accreted in the past. 

(ix) These substructures, however, may be completely 
erased if the host dwarf galaxy undergoes strong mass strip- 
ping. 

Although shells, cold stellar clumps and velocity gra- 
dients may be an indication of a past GC merger event, in 
general the detection of the associated debris will be chal- 
lenging in areas of dSphs where the field stellar population 
of the host dominates the overall light profile (see E j6,3[ ). 
Because dwarf galaxies and GCs have distinct stellar pop- 
ulations, a possible remedy may be found in future chemo- 
dynamical surveys involving hundreds of stars. These data 
sets will likely shed light on the origin of the stellar substruc- 
tures and chemodynamically distinct components detected 
in some Milky Way dSphs. 
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Table 1. Observational properties of the Fornax (Mateo 1998) and Sgr (Majewski et al. 2003) dSphs and their Globular Clusters 
(MacKey & Gilmore 2003) 



Name 


Angular sep. 


[Fe/H\ 


He 




lo gio(^) 






(kpc) 




(PC) 


(PC) 


(Lq) 


(M /p C 3) 


For dSph 


0.00 


-1.3 


400 ± 4 


2078 ± 20 


7.13 ± 0.2 


-1.14 ± 0.20 


Fl 


1.60 


-2.25 


10.0 ±0.3 


60 ±20 


4.07 ± 0.13 


0.48 ±0.07 


F2 


1.05 


-1.65 


5.8 ±0.2 


76 ± 18 


4.76 ±0.12 


1.78 ±0.07 


r o 


U.4o 


— Z.ZD 


1 ft _l_ n ft 

1.0 ± U.O 


DO ± lO 


flft 1 n m 

O.UO ± U.1Z 


Q /)7 J- n (17 

o.4 / ± U.U / 


F4 


0.24 


-1.65 


1.8 ±0.2 


44 ± 10 


4.69 ±0.24 


3.18 ±0.07 


F5 


1.43 


-2.25 


1.4 ±0.1 


50 ± 12 


4.76 ±0.20 


3.27 ± 0.07 


Sgr dSph 


0.00 


[-0.5,-1.3] 


1560 ± 20 


12600 ± 20 


7.24 ±0.2 


-2.96 ±0.20 


M54 


0.00 


-1.65 


0.91 ±0.04 


59 ±21 


5.36 ± 0.08 


4.45 ±0.05 


Terzan 7 


2.68 


-0.64 


1.63 ±0.12 


23 ±8 


3.50 ±0.10 


1.97 ± 0.07 


Tcrzan 8 


4.40 


-2.25 


9.50 ±0.72 


66 ±26 


3.67 ± 0.14 


0.72 ±0.23 


Arp 2 


3.07 


-1.65 


13.67 ± 1.85 


139 ± 49 


3.59 ±0.14 


0.35 ±0.25 
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